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Quenched crystal field disorder and magnetic liquid ground states in Tb 2 Sn 2 _:cTia ;07 
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Solid-solutions of the “soft” quantum spin ice pyrochlore magnets Tb 2 B 207 with B=Ti and Sn 
display a novel magnetic ground state in the presence of strong B-site disorder, characterized by a 
low susceptibility and strong spin fluctuations to temperatures below 0.1 K. These materials have 
been studied using ac-susceptibility and /iSR techniques to very low temperatures, and time-of- 
flight inelastic neutron scattering techniques to 1.5 K. Remarkably, neutron spectroscopy of the 
Tb^^ crystal field levels appropriate to at high B-site mixing (0.5 < a: < 1.5 in Tb 2 Sn 2 -ccTh 07 ) 
reveal that the doublet ground and first excited states present as continua in energy, while transitions 
to singlet excited states at higher energies simply interpolate between those of the end members of 
the solid solution. The resulting ground state suggests an extreme version of a random-anisotropy 
magnet, with many local moments and anisotropies, depending on the precise local configuration of 
the six B sites neighboring each magnetic Tb^^ ion. 

PACS numbers: 75.10.Dg 75.10.Kt 75.40.Gb 76.75.+i 


I. INTRODUCTION 

The rare earth pyrochlores have been a focus for the 
study of the physics of geometrical frustration, as their 
magnetic moments decorate a network of corner-sharing 
tetrahedrsPl, one of the canonical building blocks for frus¬ 
tration in three dimension^. In combination with ap¬ 
propriate magnetic couplings and anisotropies, exotic 
disordered ground states often result. In particular 
Ho2Ti2O7EJ0 and Dy2Ti2 0 7!^ exhibit net ferromagnetic 
coupling and strong local Ising anisotropy, resulting in a 
disordered spin ice ground stat^^ with Coulombic spin 
correlations and elementary excitations that behave as 
magnetic monopoleJ^^^^m. Strong quantum effects are 
not manifestly displayed by these “classical” spin ice ma¬ 
terials, and interest has now focused on several candidate 
pyrochlore magnets where a quantum analogue of spin 
ice physics may be manifest. Yb2Ti207, for example, 
possesses a spin Hamiltonian with a leading Ising spin 
component, as in classical spin ice, but with an effective 
S = 1/2 quantum moments arising from an XY-like g- 
tensor, and substantial exchange couplings that induce 
quantum dynamic^^. 

In (Ho, Dy, Yb)2Ti207, the ground state crystal field 
(GF) doublet, which describes the rare earth moment in 
its pyrochlore environment (the Ti^+ is non-magnetic), 
is well separated from the next-highest energy CF state, 
rendering all but the ground state doublet irrelevant to 


their low temperature physics. However, this is not the 
case for Tb^+ in Tb2Ti2 07, wherein the J = 6 total 
angular momentum multiplet splits into (2J + 1) = 13 
states, the lowest energy of which are a ground state 
doublet and an excited state doublet roughly 1.2 meV 
above ifP^K^. Tb2Ti207 displays antiferromagnetic in¬ 
teractions with a Curie-Weiss constant of ^ -19 
but no conventional long range order to below ^ 0.1 K, 
even though it would be expected to order magnetically 
near ^ 1 K, were its CF ground state cleanly separated 
in energy from its excited state^^^^. It displays both 
Coulomb correlations and a short range ordered spin 
ice structure at very low temperatureP^H^, as well as 
novel magnetoelastic modes t hat co uple together acous¬ 
tic phonons and CF excitation d^^ l ^^ l To explain the para¬ 
dox of Tb2Ti207, theory has invoked virtual transitions 
between the ground state CF doublet and the excited 
state doublet which stabilize a q uantu m variant of spin 
ice, precluding conventional ordeJ^^^. 

The rare-earth stanate pyrochlores, R2Sn2 07, are 
isostructural to the titanates, and have been of great 
interest, as they provide a key comparator to the frus¬ 
trated ground states of the titanates. Tb2Sn207 or¬ 
ders below Ttv ^ 0.88 K into a “soft” spin ice ordered 
phas^^. Neutron spectroscopy of Tb^+ CF states in 
Tb2Ti207 and Tb2Sn207 show the ground and lowest 
excited state wavefunctions to both be doublets sepa¬ 
rated by ^ 1.2 meV in energy. However the ground state 
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and excited state doublet eigenfunctions are made up of 
different symmetry allowed mj values in the two materi¬ 
als, primarily ^ mj=±5 in the ground state doublet and 
mj=±4 in the excited state doublefP^. The key point 
is that the Tb^+ eigenfunctions are different in these 
two sister pyrochlores, even though the low energy eigen¬ 
values are very similaili^ll^. As such, solid solutions of 
Tb 2 Sn 2 -a:Tia ,07 may constitute a novel class of spin liq¬ 
uid, emanating from disorder in the CF wavefunctions, 
which disrupt the corresponding soft quantum spin ice 
states. 

In this paper we show that strong mixing on the non¬ 
magnetic B-site in Tb 2 B 207 with B=Ti and Sn yields a 
novel disordered magnetic ground state. At variance with 
the soft quantum spin ice state of the two end members of 
the solid solution which relies on the proximity of doublet 
ground and excited states to allow virtual excitations, 
the novel disordered state at strong B-site mixing orig¬ 
inates from a disruption of the CF doublet ground and 
low lying excited states, such that these states present as 
extended continua in energy - remarkably, out to three 
times the bandwidth of these low lying CF excitations in 
either end member. It is natural to associate this with 
a distribution of local B-site environments around each 
Tb^+ ion. As the CF eigenfunctions generally determine 
the local ground state moment and anisotropy, the result¬ 
ing ground state is interpreted as an extreme version of 
a random anistropy magnet, with random ground state 
moments, on a pyrochlore lattice. 


II. EXPERIMENTAL DETAILS 

Polycrystalline samples of Tb 2 Sn 2 -a^Tia ,07 were pre¬ 
pared by standard solid state synthesis techniques. These 
same samples have be en stu died previously using ac- 
susceptibility techniqueJ^^^^. Inelastic neutron scatter¬ 
ing measurements were performed on th e SEQ UOIA 
direct geometry time-of-flight spectromete j^^ * ^^ ^ at the 
Spallation Neutron Source of Oak Ridge National Labo¬ 
ratory. Measurements were performed at T=1.5 K and 
over a wide dynamic range of energies, using incident 
neutron energies, E^, of 11 meV, 45 meV and 120 meV, 
in order to probe the appropriate CF excitations and low 
energy spin dynamics. 

/iSR measurements were performed at the Paul Scher- 
rer Institute using the GPS spectrometer for tempera¬ 
tures greater than 1.6 K and the LTF spectrometer for 
low temperatures below 1.6 K. The muon decay asym¬ 
metry was recorded as a function of time in all samples 
in a small longitudinal field (50 Oe) in order to decouple 
possible low field muon sites. 


III. NEUTRON SCATTERING RESULTS 

Figure shows colour contour maps of the low energy, 
inelastic scattering for six different Tb 2 Sn 2 -a:Tia ,07 sam- 
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FIG. 1. (color online) The measured S(|Q|,/ia;) for 

Tb 2 Sn 2 -ccTb 07 at T=1.5 K is shown. These low energy, 
high resolution data was taken with Ei=ll meV neutrons. All 
data have been corrected for empty can background and de¬ 
tector efficiency. Well defined transitions between the ground 
state crystal field doublet and the lowest excited crystal field 
doublet near 1.5 meV are observed for the end members 
Tb 2 Sn 207 , X = 0, and Tb 2 Ti 207 , x — 2^ while intermedi¬ 
ate concentrations see a much broader bandwidth to these 
low energy excitations. 


pies from x = 0 to x = 2, using E^=ll meV incident 
neutrons. Well defined low energy CE excitations for the 
end members, x = 0 and x = 2, appear as sharp peaks in 
energy between 1-2 meV with weak dispersion, such that 
they display a minimum near \Q\ ^ 1. 25 A ^ for both 
case ji3 | i5H iIj^ However, as seen in Eig. [l] the well defined 
CE excitations are lost at all wavevectors for intermedi¬ 
ate concentrations, with some enhanced spectral weight 
concentrated in the 1-2 meV range for the most lightly 
disordered x = 0.2 sample. Remarkably for x = 0.5, 1 
and 1.5, we observe a very broad distribution of low en¬ 
ergy magnetic scattering at all \Q\s studied extending to 
three times the energy of the lowest lying CE excitations 
in either end member. 

Eigure shows cuts of these and related data sets, 
taken with different incident energies (E^s) and over dif¬ 
ferent energy ranges, to focus on different CE excitations 
between 1 meV and 50 meV. Eigure (a) shows inten¬ 
sity vs energy cuts for the E^=ll meV data of Eig. 

It shows the excitations from the ground state doublet 
to the excited CE doublets near 1.5 meV in the two end 
members, as well as the absence of these transitions in 
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FIG. 2. (Color online) Cuts of the inelastic neutron scatter¬ 
ing data for Tb 2 Sn 2 -ccTiic 07 samples at T=1.5 K are shown, 
(a) and (b) show inelastic scattering taken with Ei=ll meV 
and 45 meV neutrons, respectively, (c) shows inelastic scat¬ 
tering taken with £^=120 meV neutrons. All data have been 
corrected for empty can background and detector efficiency. 
The high energy crystal held excitations in (c) were ht to 
extract their energies, and these are plotted in the inset to 
(c). Data sets in (a) and (b) have been integrated between 1 
< \Q\ < 1.4 A“^, while that in (c) has been integrated 
between 2 and 3 A' ^ in |Q|. 


the solid solutions. The x = 0.5, 1.0, and 1.5 samples all 
show much enhanced quasi-elastic scattering at low en¬ 
ergies (<0.8 meV) compared with the end members, and 
the spectral weight of the diffuse scattering in the regime 
of the lowest lying CF excitation is now extended to ^ 6 
meV, consistent with the colour contour maps of Fig. 

At higher energies, i.e. the E^=45 meV data in Fig. 

(b) and E^=120 meV data in Fig. (c), we observe a 
continuous evolution of the CF excitations with concen¬ 
tration. In the two parent materials the CF excitations 
in Fig. (b) have recently been shown to correspond 
primarily to transitions from the doublet ground state to 
singlet CF states at ^ 10, 14 and 17 meV for Tb 2 Ti 207 , 
and ^ 10 and 15.5 meV for Tb 2 Sn 2 07 , while those in 
Fig. i (c) correspond primarily to transitions from the 
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FIG. 3. (Color online) The real, and imaginary, y// parts 
of the ac susceptibility are shown for Tb 2 Sn 2 -ccTia ,07 at low 
temperatures < 1 K. Measurements were performed in zero 
DC magnetic held, at a frequency of 500 Hz. 


doublet ground state to singlet CF states at ^ 49 meV for 
TbaTiaOy, and ~ 34 and 34.8 meV for TbaSnaOT^i^. At 
these higher energies, where transitions to singlet excited 
states are primarily studied, we observe a remarkably 
smooth progression in the B-site disordered systems, as 
the appropriate CF excitations evolve continuously from 
the X = 0 to X = 2 limit, as seen from ^ 32 meV for 
Tb 2 Sn 207 to ^ 49 meV for Tb 2 Ti 207 in Fig. [^(c). The 
excitations were fit to Lorentzians to extract the CF ex¬ 
citation energy as a function of concentration, and this 
is shown in the inset to Fig. [^(c). Clearly the evolution 
of this high energy CF excitation is linear with compo¬ 
sition. Furthermore, while the CF excitation broadens 
somewhat for intermediate compositions, it remains well 
defined and at energies which are different from that of 
the end members, very distinct from the transitions be¬ 
tween the doublet ground state and doublet first excited 
states shown in Fig. [^and Fig. |^(a). 


IV. AC SUSCEPTIBILITY AND fiSK 
MEASUREMENTS 

At lower temperatures, the phase behavior and low 
energy spin dynamics across this series of solid solutions 
were studied using ac susceptibility and /iSR techniques. 
The real, %/, and imaginary, %//, parts of the susceptibil¬ 
ity, measured in zero external field and at / = 500 Hz, 
are shown in Fig. |^(a) and (b), respectively, for temper- 
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FIG. 4. (Color online) The /xSR relaxation rate measured 
with a small 50 Oe longitudinal field is shown as a function of 
temperature for selected Tb 2 Sn 2 -ccTR 07 samples. Both the 
relaxation rate and the temperature are plotted on logarith¬ 
mic scales. 


atures below 1 K (some of these data and methods have 
been previously published in Ref l3Q)) . For Tb 2 Sn 2 07 we 
confirm the peak in y/ and concomitant inflection in y// 
at Tat ^ 0.83 K, signifying the phase transition to the 
“soft” spin ice ordered phase. These features are wiped 
out for all subsequent levels of Sn/Ti mixing, indicat¬ 
ing the fragility of the ordered phase in Tb 2 Sn 2 07 l^. 
At the other end of the solid solution, x = 2, we ob¬ 
serve a weak peak in both y/ and y// indicating glassy 
behavior for temperatures of ^ 0.23 K and below, con¬ 
sistent with earlier studie^^^K^. These features display 
little frequency dependence in contrast to expectations of 
a canonical spin glass transitioiP^. However, for strong 
B-site disorder, in particular x = 0.5, 1, and 1.5, only 
very broad features with low absolute values for either y/ 
or y//, and no tendency of y/ to go to zero as temperature 
approaches zero, are observed. This is again consistent 
with a highly disordered liquid-like ground state distinct 
from either end member. 

/iSR measurements in a small longitudinal field (50 Oe) 
were performed on this same solid solution series, span¬ 
ning the temperature range from 20 mK to 300 K. We 
observed a nearly single exponential relaxation to the de¬ 
cay asymmetry over the entire temperature regime and 
for all samples. The resulting relaxation times are shown 
as a function of temperature in Fig. A weak inflec¬ 
tion is observed near the transition to the “soft” spin ice 
ordered phase in Tb 2 Sn 207 , in t he temperature depen¬ 
dence to its relaxation rate^^®. Both Tb 2 Sn 207 and 
Tb 2 Ti 207 display a high temperature plateau between 
^ 10 K and 200 K, due to the thermal occupation of the 
lowest lying CF excited states, ^ 15 K above the ground 
state for either. The relaxation rate in Tb 2 Ti 2 07 grad¬ 
ually rises to ^ 2.5 /isec“^ at the lowest temperatures 
measured. 

The /iSR signature for spin freezing in the “ordered 



FIG. 5. (Golor online) The low temperature muon relax¬ 
ation rates for different members of the Tb 2 Sn 2 -xTia ,07 se¬ 
ries. Note that both axes are logarithmic. The inset shows dif¬ 
ferences of the relaxation rates measured in a 50 Oe field, on a 
linear scale, following protocols of a 800 Oe field-cooling (FG) 
and a zero field cooling (ZFG) for samples with x = 0.1, 0.2 
and 1.0, compared with x = 0 (from Ref. [40l) . As with the 
main panel, the temperature axis is logarithmic. Lines are 
guides to the eye. 


spin ice” ground state in Tb 2 Sn 207 proved to be quite 
subtle. None of the usual signs of spin freezing, such as 
recovery of a 1 /3rd tail of the asymmetry or strong sensi¬ 
tivity to an appli ed lon gitudinal field, could be detected 
below Tc ^ 0.9 KpSESl, Indeed the clearest indication for 
frozen spin correlations was found to be the difference of 
the relaxation rates measured under field-cooling (FC) or 
zero field cooling (ZFC) protocols. We carried out these 
same FC/ZFC protocols for the Tb 2 Sn 2 -a:Tia ,07 series 
with X = 0.1, 0.2 and 1, in addition to x = 0. As shown 
in Fig. this difference, and therefore the ordered spin 
ice state, is rapidly suppressed at finite x (see inset of 
Fig. [^. Both the amplitude of the FC/ZFC difference 
in relaxation rate, as well as its onset temperature, are 
reduced from 1.2 K at x = 0 to ^ 0.5 K for x = 0.1, a 
temperature which corresponds to a broad maximum in 
the ac susceptibility (Fig. [^(a)). This FC/ZFC signature 
completely disappears for x = 0.2 and x = 1.0. While y' 
displays a peak as a function of temperature at ^ 0.25 K, 
possibly indicative of a spin glass transition, only weak 
features are observed for 0.2 < x < 1.8. We therefore 
conclude that the ground state for Tb 2 Sn 2 -ccTia ,07 and 
0.2<x<1.8isa fully dynamical spin liquid. 

For intermediate concentrations, the high temperature 
plateau fills in (see inset to Fig. |^, consistent with the ef¬ 
fect of strong B-site disorder on the low lying crystal field 
levels, shown in Figs and Again at low temperatures 
the relaxation rate rises gradually to a low temperature 
limit of ^ 2 /isec“^, showing a progressively smaller low 
temperature limit than either end member; again consis¬ 
tent with a disordered and dynamic ground state, and 
no hysteretic behavior, in the presence of strong B-site 
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FIG. 6. (Color online) A minimal model for the effect of 
B-site disorder on the low energy CF levels in Tb 2 B 2 07 is 
illustrated. B-site disorder can split the doublet CF states, 
as Tb3+ 

is a non-Kramers ion, but the singlets are simply 
shifted in energy. Disordered B-site samples possess a distri¬ 
bution of the seven variants (3 Sn^^, 3 Ti^^ ions, 4 Sn^^, 2 
Ti^^ ions, etc.) of local CF environment at the Tb^^ site. 
Transitions from the ground state to excited CF states are 
indicated with red arrows, while those terminating with a red 
dot indicate inelastic peak positions expected for high B-site 
disorder, due to averaging over local Tb^^ environments. The 
local environment at the Tb^^ site in the pyrochlore structure 
is indicated at the top. For clarity, the six ions around 
each Tb^^ ion are omitted. 


disorder. 


V. DISCUSSION AND CONCLUSIONS 

The Tb3+ ion possesses an even number of electrons, 
in contrast to Dy^+ and Yb^+, and thus its CF config¬ 
uration is not protected by Kramers theorem. The 13 
CF levels in Tb 2 Ti 207 and Tb 2 Sn 2 07 then appear as 
both singlets and doubletJ^. Remarkably, our neutron 
spectroscopy of the CF states in the presence of B-site 
disorder strongly suggests that doublet states present as 
extended continua in energy, while singlet states have 
their energy eigenvalues and lifetimes modified, but ap¬ 
pear less strongly affected. This is consistent with a min¬ 
imal model of the low lying CF levels, where CF doublets 
split in the presence of B-site disorder, but singlets sim¬ 
ply have their eigenvalues modified, as illustrated in Fig. 
^ As the local crystalline environment around a Tb^+ 


ion, shown at the top of Fig. is determined by a ring 
of B-site ions and associated ions, a strongly dis¬ 
ordered Tb 2 Sn 2 -a^Tia ,07 sample would display the ap¬ 
propriate statistical distribution of environments for the 
Tb^+ sites, corresponding to all permutations of Ti^+ 
and Sn^+ ions decorating the six neighboring B-sites. 
The B-site ions have different ionic radii, and would also 
perturb the local positions slightly. The CF eigen¬ 
functions and eigenvalues determine both the local Tb^+ 
ground state moment and its anisotropy. Such a distri¬ 
bution of quenched, disordered CF environments would 
give an extreme version of a random anisotropy mag¬ 
net, in the presence of random moment sizes, which, at 
a minimum, would also disorder the dipolar part of the 
spin Hamiltonian in Tb 2 Sn 2 -ccTia, 07 . A related effect on 
the doublet CF ground state of Pr^+ in A-site disordered 
Pr 2 -a:Bia,Ru 207 has likely been observed^. 

Taken together, these experimental results show a re¬ 
markable route to a new type of strongly fluctuating mag¬ 
netic liquid state, based on strong quenched disorder of 
the Tb^+ CF environment. Given that both the Sn'^^ 
and Ti^+ sites are non-magnetic, one may naively expect 
the mixing of these elements on the B-sites of Tb 2 B 207 
to have little effect on the low temperature physics dis¬ 
played across the solid solution. Yet it has a dramatic 
effect; it wipes out the subtle transition to a “soft” spin 
ice ordered phase in Tb 2 Sn 207 , and produces a fluctuat¬ 
ing spin liquid state at intermediate concentrations. 

These results also underline the crucial role played 
by the low lying excited CF doublet for Tb^+ in both 
end members, x = 0 and x = 2. The modification of 
the ground state and lowest excited state doublet wave- 
functions on moving from x = 0 to x = 2, implies a 
disruptive change to the virtual excitations between the 
ground state and low lying excited states. As this is the 
mechanism proposecP^ by which quantum fluctuations 
are introduced to Tb 2 Ti 2 07 , and which presumably are 
also at play in Tb 2 Sn 2 07 , the advanced characterization 
of the Tb 2 Sn 2 -a:Tia ,07 solid solutions provides an im¬ 
portant perspective on such quantum fluctuations in the 
presence of relevant disorder. 
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